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We investigate the behavior of the magnetic pressure, , in fully turbulent MHD 
flows in "1+2/3" dimensions by means of its effect on the probability density 
function (PDF) of the density field. We start by reviewing our previous results for 
general polytropic flows, according to which the value of the polytropic exponent 
determines the functional shape of the PDF. A lognormal density PDF appears in 
the isothermal (7 = 1) case, but a power-law tail at either large or small densities 
appears for large Mach numbers when 7 > 1 and 7 < 1, respectively. In the 
isothermal magnetic case, the relevant parameter is the fleld fluctuation amplitude, 
5B/B. A lognormal PDF still appears for small field fiuctuations (generally the 
case for large mean fields), but a significant low-density excess appears at large 
fluctuation amplitudes (weak mean fields), similar to the behavior at 7 > 1 of 
polytropic flows. We interpret these results in terms of simple nonlinear MHD 
waves, for which the magnetic pressure behaves linearly with the density in the 
case of the slow mode, and quadratically in the case of the fast wave. Finally, we 
discuss some implications of these results, in particular the fact that the effect of 
the magnetic field in modifying the PDF is strongest when the mean field is weak. 



1 Introduction 

A fundamental feature of compressible turbulence is the formation of den- 
sity fluctuations, a property of central interest in astrophysics, as density 
excesses ("clouds") in the turbulent interstellar mfidium (ISM) exhibit nu- 
merous statistical properties over a range of sizescrEI whose origin is not yet 
well understood. Although a full understanding of star formation requires 
knowledge of the full (multiple-point) statistics in order to determine mean 
densities as a function of region size, a first step towards this goal is the de- 
scription and physical understanding of the one-point statistics, or probability 
density function ("PDF") of the mass density field. 

Previous studies have shown that the density PDF in_turhiilent compress- 
ible flows has a lognormal form in the isothermal casejOO'EfEl but exhibits a 
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power-law tail at densities larger (smaller) thaii.the.mean for flows with poly- 
tropic exponents smaller (larger) than unityOO'El It should_be_poted that 
some of those works referred to purely hydrodynamic flows jlfl'O while the 
rest referred to magnetOj- hydrodynamic (MHD) flows. 

In particular, ref. p3] presented a heuristic model for the development of 
the lognormal and power-law PDFs, based essentially on the behavior of the 
speed of sound with density in those types of flows. In this sense, the density 
PDF is a diagnostic for the dependence of the pressure with density. However, 
such model applied only to non-magnetic flows, while the real ISM is most 
likely magnetized to a signiflcant extent. In this paper we present prelimi- 
nary results on the density PDF of turbulent magnetized flows in "1-1-2/3" 
dimensions, as a first attempt to characterize the behavior of magnetic pres- 
sure with density in the fully turbulent regime. Previous works have focused 
on the pressure produced by weakly nonlinear Alfven wavesOilj, but here we 
consider fully turbulent regimes with arbitrarily large magnetic fluctuation 
amplitudes. 

2 The Non-Magnetic Case 

The form of the density PDF in non-magnetic flows and its relation to the 
effective equation of state of the system has been_understood in terms of a 
heuristic model by Passot and Vazquez-Semadenijid herafter PVS98 (see also 
ref. [B]). As in ref. this model idealizes the generation of turbulent den- 
sity fluctuations as a "multi-step" process, in which the local density at any 
given point in the flow is the result of a series of jumps due to the continuous 
passage of shock waves. The generation of densities is an iterative multiplica- 
tive process, in which every new density is obtained through a jump from the 
previous one, giving for the final density pt in terms of the intermediate steps 
Pi'- Pi = PoY\i{pi+i/ Pi) ■ In the isothermal case, the density jump is given 
by M^, where M is the Mach number of the shock and, because the speed of 
sound is constant, a given velocity jump is always characterized by the same 
Mach number, regardless of the local density, so that individual jumps can be 
regarded as independent, but extracted from the same distribution. In terms 
of the variable s = In p, the iterative process is additive so that the Central 
Limit Theorem can be applied in the limit of a large number of jumps, lead- 
ing to a normal distribution for s and thus a lognormal distribution for p, as 
observed in the numerical experiments mentioned above. 

Concerning the variance of the distribution, PVS98 have suggested (and 
confirmed numerically) , through an analysis of the shock and expansion waves 
in the system, that for a large range of Mach numbers the typical size of the 
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logarithmic jump is expected to be o"s ^ Mrmg, where A/i., 
number. The mean sq of the distribution can be directly evaluated from the 
mass conservation condition (p) = e^P{s)ds = 1, where P{s) is the PDF 
of s, yielding sq = —a'^/2. The isothermal model PDF for s thus reads 

with (T^ = l3M^^g and /? a proportionality constant of order unity. 

In the general poly tropic case where the pressure P behaves as P oc p''', 
with 7 the polytropic exponent, the local Mach number M(s) at a density 
p = is related to the one at the mean density (M) by M{s) = M exp'^'^"^)*/^, 
suggesting an ansatz where the PDF keeps the same dependence on M, pro- 
vided the above replacement is made. After relocating the term in so from 
inside the exponential function to the normalization constant, the model PDF 
for the polytropic case reads 

r_s2g(7-l)s -, 

P{s; 7)ds = C(7) exp |^ — 01(7)3] ds. (2) 

This equation shows that when (7 — l)s < 0, the PDF asymptotically ap- 
proaches a power law, while in the opposite case it decays faster than a lognor- 
mal. Thus, for < 7 < 1, the PDF approaches a power law at large densities 
(s > 0), and at low densities (s < 0) for 7 > 1. 

3 The Magnetic Case 

In what follows we restrict ourselves to the isothermal case (7=1) and to a 
propagation along a uniform ambient magnetic field Bq and concentrate on the 
deviations from the corresponding lognormal PDF induced by the magnetic 
field. 



3.1 Numerical Results 

The numerical simulations solve the MHD equations using a pseudo-spectral 
method in a slab geometry ("1+2/3D") (variabihty is considered only with 
respect to the spatial variable x for the three components of the velocity and 
magnetic fields). A resolution of 2048 grid points is used allowing to handle 
large enough Mach and Reynolds numbers with only regular second-order 
viscosity. A random acceleration is applied on the y- and ^-components of the 
velocity field, i.e. transverally to Bo, thus inducing Alfven waves into the flow. 
Since these waves have finite amplitudes, they in turn induce fast and slow 
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Figure 1. Density PDFs for two ID simulations of polytropic turbulence, one with 7 = 0.3 
(left), and the other with 7 = 1.7 (right), both with M = 3. 



magnetosonic waves. The simulations are evolved over long times (several 
tens to a few hundred crossing times at the rms flow velocity) in order to 
obtain meaningful statistics for the density PDFs. 

The simulations are essentially characterized by the sonic and Alfvenic 
Mach numbers, respectively defined as Mg = u/cg and Afa = m/wa, where u 
is the rms velocity, Cg the sound speed and va = B^/ p the Alfven speed. In 
order to investigate the effect of the magnetic field exclusively, we consider two 
runs (denoted I and II) with approximatively the same value of Afg (4.06 and 
3.74 respectively) but very different values of A/a (0.36 and 1.66 resp.). These 
quantities are integrated over the duration of the run, excluding the first few 
time units in order to avoid including the uniform-density initial conditions. 
The rms relative field fiuctuations 5B/B are for runs I and II respectively 
0.75 and 6.65. 

In fig. H we show the PDFs for the two magnetic runs, together with 
lognormal fits (dashed lines). It is clearly seen that, contrary to what one 
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Figure 2. Density PDFs of two magnetic simulations with similar values of the sonic Mach 
number Ms but different values of the Alfvenic Mach number Ma. The simulation with 
small Ma. (large mean field) is seen to have a nearly lognormal PDF, while the one with 
large Ma (weak mean field) is seen to have a large excess at low densities, indicative of an 
effective "magnetic" polytropic exponent larger than unity. 



might expect, the high-Ma run, which should be closest to the pure-hydro 
case, exhibits a strong deviation from lognormality, while the low-A/a PDF is 
quite close to a lognormal. This implies that the case with the weakest mean 
field is the one with the strongest effect of magnetic pressure on the density 
PDF. 

3.2 Interpretation in Terms of Simple MHD waves 

A preliminary interpretation of the behavior reported in the previous sectipn 
can be given in terms of the so-called "simple MHD waves" (see, e.g., ref. p]). 
These are the finite-amplitude equivalent of linear MHD waves, and have the 
same three well-known modes: Alven, slow and fast. A simple wave refers 
to a solution that depends only on a single variable, any combinations of the 
spatial (x) and temporal (t) independent variables. However, as mentioned 
above, simple waves have finite amplitudes, and in general they develop shocks 
in a finite time. 

A number of properties ^f simple waves has been reported by Manno 
(see also Jeffrey and Taniuticl). Among them, most relevant for the present 
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discussion are the propagation velocities of the three modes, and the relations 
between the density and magnetic field fluctuations. The Alfven mode with 
speed V such that = va^ is not associated with density fluctuations. The 
fast and slow speed are given by 



Vl= : \l±\l- 



r} (3) 



and the relation between the magnitude of the magnetic field b and the density 
is given by 

dp b ' ^ ' 

where va^ = b^/p, = b'^/p^ a-nd Cg is the sound speed. An important point 
to note is that eqs. (||) and imply a positive correlation between fa-and p 
for the fast mode (V+) and an anticorrelation for the slow mode (VL),a. 

Inserting eq. (^) in eq. one can find the density dependence of the 
magnetic pressure. Excluding the case where the /3 of the plasma is of order 
unity with at the same time small to moderate field distorsions, this depen- 
dence siinplifies and reads (general numerical solutions have been presented 
by MannH): 

&2 c^p 

-2-wai-a2-^ (slow mode) (5) 
b^ K. 2a3p2 (fast mode), (6) 

where b^ = b'^ + b\, with b^ = by + b^ being the magnitude of the magnetic 
field fluctuation. Note that in our 1+2/3 geometry, bx{= Bo) is constant. The 
quantities a; denote integration constants. 

From these relations, it can be seen that the magnetic pressure, oc b^, 
scales roughly linearly with the density (albeit inversely) in the case of the 
slow mode, but quadratically in the case of the fast mode. Note also that, in 
the case of small field fiuctuations (most often the case if /? = Cg/(6^/p) <C 1) 
the coefficient of p in eq. (^) is small, implying that, for the slow mode, 
large density fluctuations can occur even for small variations of b. Thus, the 
slow mode is expected to dominate density fluctuation production in the case 
of small field fluctuations. Instead, at large enough field fluctuations, the 
quadratic dependence of b^ on p of the fast mode eventually overwhelms the 
linear dependence of the slow mode (and also of the thermal pressure), so the 
fast mode is expected to dominate the density fluctuation production. 

The dependence of the magnetic pressure on the density seems to describe 



the observations of sec. |3.1| adequately because, according to the results of the 
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non-magnetic case, a pressure that depends linearly on the density produces 
a lognormal density PDF, while one with a higher-than- linear dependence 
produces a near power law at low densities. 

4 Conclusions 

Summary 

In this paper we have reviewed our previous resultsS on the development 
of the density PDF as a consequence of the effective (polytropic) equation of 
state, and applied them to an understanding of the. PDF in the magnetic case. 

In the context of the model presented in ref . [U] , isothermal non-magnetic 
flows have lognormal mass density PDFs while in the non-isothermal cases, 
a power-law tail develops at high Mach number, at either p > (p) or p < (p) 
depending on whether 7 < 1 or 7 > 1, respectively. 

In the magnetic case, we considered only isothermal cases and a propa- 
gation parallel to the ambient magnetic field, but reported that a deviation 
from the lognormal PDF occurs when the magnetic fluctuations are large, a 
case expected when the mean field is small. We interpreted this effect in 
terms of so-called "simple" , finite-amplitude nonlinear MHD waves. Indeed, 
for the slow mode of these waves, the magnetic pressure behaves linearly 
with the density p, while for the fast mode, the magnetic pressure behaves 
quadratically with p. Together with the fact that for the slow wave the den- 
sity is weighted by a small factor, this suggests that the production of density 
fluctuations is dominated by the slow waves at small field deviations (i.e., 
Bo large) giving a lognormal PDF again, while for weak mean fields (large 
field deviations) the quadratic density dependence of &^ produces an excess 
at small densities in the PDF, corresponding to the 7 > 1 case of the poly- 
topic description. A more detailed study, including in particular the case of 
non-parallel propagation is in progress. 

4-2 Astrophysical Implications 

The results of this paper have a number of implications in the astrophysical 
context. First, we have seen that the main parameter determining the ef- 
fect of the magnetic field appears to be the magnetic fiuctuation amplitude, 
6B/B, rather than (3 alone, which is the parameter most frequently used to 
characterize MHD flows. In other words, it is important to have a knowl- 
edge of the relative importance of the turbulent fluctuations (as given by the 
Alfvenic Mach number) in addition to simply the ratio of thermal to magnetic 
pressures (as given by /3). 
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Second, as far as the density PDF is concerned, the effect of the magnetic 
field is most notorious when the field fluctuations are large (generally, when 
the mean field is weak) . This suggests that the limit of vanishing field does 
not approach the non-magnetic case, and in this sense the latter is singular. 
This may imply that it is inadequate to model the ISM as a non-magnetic flow 
even if the degree of magnetization is very low, as has staiJted to be recently 
contended at the level of molecular clouds (see, e.g., ref. [c3]). 
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